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ABSTRACT 



We present the first detection of parallax effects in a gravitational microlensing event. 
Parallax in a gravitational microlensing event observed only from the Earth appears as a 
distortion of the (otherwise symmetrical) lightcurve due to the motion of the Earth around 
the Sun. This distortion can be detected if the event duration is not much less than a year 
and if the projected velocity of the lens is not much larger than the orbital velocity of 
the Earth about the Sun. The event presented here has a duration (or Einstein diameter 
crossing time) of t = 220 days and clearly shows the distortion due to the Earth's motion. 
We find that the projected velocity of the lens is v = 75 ± 5km/s at an angle of 28° ± 4 
from the direction of increasing galactic longitude, as expected for a lens in the galactic 
disk. 

A likelihood analysis of this event yields estimates of the distance to and mass of the 
lens: D\ ens = 1.7^Q^kpc and M = 1.3^QgM Q . This suggests that the lens is a remnant 
such as a white dwarf or neutron star. It is possible, though less likely, that the lens is a 
main sequence star. If so, we can add our upper limit on the observed flux from the lens to 
the analysis. This modifies the estimates to: D\ ens = 2.8*Qg kpc, and Mi ens = 0.6^Q2M Q . 
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Introduction 



In less than two years since the first candidate gravitational microlensing events were 
discovered (Alcock et al. 1993, Aubourg et al. 1993, Udalski et al. 1993), microlensing 
has been confirmed spectroscopically (Benetti et al. 1995) and more than 60 microlensing 
events have now been discovered (Udalski et al. 1994, Alcock et al. 1995a, 1995b, Bennett 
et al. 1995). The vast majority of the events have been discovered towards the Galactic 
bulge, and the total microlensing optical depth seen towards the bulge seems to be a factor 
of ~ 3 larger than predicted (Griest et al. 1992, Paczyhski 1992). These results suggest 
that standard Galactic models are probably in need of revision (Paczyriski et al. 1994, Zhou 
et al. 1994, Gould 1994), and this may have important implications for the interpretation 
of microlensing results toward the LMC (Alcock et al. 1995a, Aubourg et al. 1995). 

Microlensing light curves are generally distinguished from other types of stellar vari- 
ability by their very simple form. This is a great benefit when trying to detect these 
very rare events, but it makes their interpretation more difficult because a wide variety of 
lensing events can give rise to very similar light curves. In this paper, we present the first 
detection of the microlensing parallax effect which, as we shall see, offers great promise for 
resolving the degeneracies in microlensing events. 

In most microlensing events, it is a very good approximation to assume that the source 
star is point-like, the lens mass is also point-like, and that the velocities of the source, the 
lens, and the observer are all constant in time. Under these assumptions, microlensing 
light curves take a very simple form: the amplification is given by 

A(t) = -^±1=; u(t) = ^4 + [2(t - t )/m (1) 
MV u +4 v 

where t is the Einstein diameter crossing time defined by 

i/2 = R E /v t ; R\ = 4GMx(l - x)L/c 2 , (2) 

where Re is the Einstein ring radius, vt is the transverse velocity of the lens relative to 
the (moving) observer-source line, M is the lens mass, L is the observer-source distance, x 
is the ratio of the observer-lens and observer-source distances, and uq = b/RE where b is 
the closest approach of the lens to the observer-source line. We can see from eqs. (1) and 
(2) that A[t) is described by only three parameters: uq, to, and t (in practice, there is an 
additional parameter for the unlensed flux of the star in each passband). Furthermore, the 
parameters uq and to reveal only the relatively uninteresting information about when the 
lens came closest to the line of sight and how close it came. The information about M, x, 
and vt, which would help us determine whether the lensing object is part of a bulge or disk 
population, is all folded into the single parameter t, so it is impossible to determine these 
parameters separately (though probability distributions may be obtained, e.g. Griest et 
a/.1992). 
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The situation is improved for events in which one or more of the assumptions leading 
to the generic 3-parameter light curves described in eqs. (1) and (2) are violated, (e.g., 
Mao & Paczyhski 1991; Gould 1992, 1994; Gould et al. 1995). The MACHO collaboration 
has observed two such events in the first year bulge data. One is the event presented in 
this paper and the other is the binary lensing event first seen by OGLE (Udalski et al. 
1994b). Our data on the binary event resolves the finite size of the source star (Bennett et 
al. 1995), and when our analysis of this event is complete, it should yield two constraints 
on combinations of M, x, and the lens velocity. 

The situation is similar for the event described here, which is the longest of 45 events 
in our first-year bulge data. The light curve for this event is shown in Fig. 1. The dashed 
curve is the best-fit symmetric light curve described by eqs. (1) and (2). Clearly, the 
light curve of this star exhibits a significant asymmetry and has a large deviation from the 
best-fit symmetric light curve. However, the light curve is still achromatic, so it is likely 
to be a microlensing event. Such an asymmetry could be caused by a deviation from the 
constant velocity assumption for the source, lens, or observer. In principle, it is impossible 
to strictly distinguish between these three possibilities because each can contribute equally 
to vj, which is the transverse velocity of the lens with respect to the line of sight to the 
source star. However, since the orbit of the Earth is known, one can attempt to fit the light 
curve under the assumption that the only significant deviation from constant velocities is 
the orbit of the Earth. If a good fit is obtained assuming the known period, orientation 
and phase of the Earth's orbit, then it is reasonable to assume that the orbit of the Earth is 
indeed responsible for the deviation of the light curve because it is quite unlikely that the 
orbital parameters of the source star or the lens would exactly match those of the Earth. 
As we shall see, the event presented in this paper is a case where the orbital motion of 
the Earth is probably responsible for the deviation of the lightcurve. This is a fortunate 
circumstance as it allows us to compare the projected Einstein ring diameter crossing time 
with the known size of the Earth's orbit and thereby obtain a second constraint on the 
three unknown parameters of the lens, M, x, and vt- 

The star shown in Fig. 1 is located at a = 18 h 03 m 34.05 s , 5 = -28°00 / 18.9" (J2000) 
which in ecliptic coordinates is A = 270.85° and j3 = —4.65°. The small value of is a 
potential difficulty for our attempt to solve for the mass and location of the lens because 
the sign of the perpendicular component of the projected transverse velocity becomes 
undetermined in the (3 = limit. To include the orbital motion of the Earth in eq. (1), 
we must replace the expression for u(t) with 

u 2 (t) =Uq + co 2 (t- t ) 2 + a 2 sin 2 [fi(£ - t c )} 

+2asin[fi(£ - t c )] [u(t - t ) sin (9 + u cos6] (3) 
+a 2 sin 2 (5 cos 2 [fi(£ — t c )] + 2a sin (5 cos[fi(£ — t c )] [u(t — to) cos 6 — uq sin 6*] 

where 9 is the angle between and the North ecliptic axis, u = 2/t, and t c is the time when 
the Earth is closest to the Sun-source line. Note that uq is no longer the minimum distance 
between the lens and the observer-source axis as it was in the constant velocity case, but 



4 



is the minimum distance between the lens and the Sun-source line. The parameters a and 
O are given by 

a = w(lAU) _ ecQs[Qo{t _ )]} (4) 

and 

ft(f - t c ) = n (t- t c ) + 2e sin[fi (^ - t p )] , (5) 



where t p is the time of perihelion, v = vt/(l — x) is the transverse speed of the Macho 
projected to the position of the Sun, fio = 27ryr _1 , and e = 0.017 is the eccentricity of the 
Earth's orbit. 

The dashed curve in Fig. 1 is the best-fit constant velocity microlensing light curve, 
and the solid curve is the fit including the motion of the Earth. The fit parameters and x 2 
values are listed in Table 1. As seen from Table 1, the x 2 P er degree of freedom is reduced 
dramatically from 8.42 to 1.54 when the terrestrial orbit terms are included. Formally, a 
X 2 per d.o.f. of 1.54 does not indicate a good fit for 206 degrees of freedom, but this is 
not unusual in our data for stars of this magnitude. If the 2 red band outlier points at 
days 135 and 201 are (arbitrarily) removed from the fit, the x 2 P er d.o.f. drops to 1.08. In 
fact, these two outlier measurements do have unusually large point-spread- function fit x 2 
values suggesting suspect photometry, but we do not impose an automatic cut on the PSF 
X 2 value for stars that are as bright as this. 



fit# 


to 
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jj (kffi) 

V sec > 
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x 2 
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143.7(8) 
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1.000(3) 


1.000(4) 
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129.0(6) 
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0.101(5) 


9.93 


1.012(3) 


1.001(4) 


oo 


0.0 


1752 



Table 1: Fit results for 213 measurements in both passbands. Fit # 1 is the full fit 
including parallax effects. Fit # 2 is the usual "constant velocity" fit, which is sufficient 
to describe most microlensing light curves. The error estimates given are the maximum 
extent of the surface in parameter space, which has x 2 greater than the best-fit value by 
1. 



The direction of the best-fit velocity, 8 = —1.01 rad from ecliptic North, is 28° away 
from the direction of increasing galactic longitude. This is what we expect for a lens 
located in the galactic disk, where much of the velocity should be due to the disk rotation, 
and thus the lens is 'overtaking' the line connecting the moving Sun to the galactic bulge 
which is stationary (on average). 

Perhaps the most interesting aspect of this microlensing parallax observation is that it 
allows us to learn much more about this event than we can learn about 
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ordinary microlensing events. Substituting vt = v(l — x) into eq. (2), we obtain 



M(x) 



I - x v 2 P c 2 



(6) 



x 16GL 



Fig. 2 shows M(x) for the values given in Table 1 (fit # 1) assuming that L = 8kpc. 
We can see from Fig. 2 that the lens could be either a brown dwarf in the Galactic bulge, 
an M dwarf between 2 and 6 kpc distant, or a more massive star nearby. If we assume 
that the star is currently burning its nuclear fuel, then we can also constrain this latter 
possibility because such a star would contribute additional unlensed flux to the observed 
light curve, since it would be within a few milli-arcseconds of the lensed star. 

We have redone our parallax event fit including the possibility that there is some 
unlensed light superimposed upon the lensed star, an R m 15.9 clump giant. The best fit 
brightness of this assumed companion is — 4±4% of the normal brightness of the lensed star. 
We will take ± 4% as our upper limit on the apparent brightness of the lens relative to 
the source. Then, assuming an .R-band mass-luminosity relation of the form L(M) = M 2 6 
(Mihalas & Binney, 1981) in the mass range 0.6M Q < M < M Q , and allowing for 1.5 
magnitudes of extinction, we obtain an upper limit on the mass of a main-sequence lens 
star, which we will use in Fig. 3 below. 

Another, somewhat more general constraint on x and M can be obtained if we make 
use of our knowledge of the velocity distributions of the source and lensing objects, since 
the likelihood of obtaining the observed value of v is a strong function of the distance to 
the lens. For example, if the disk and bulge velocity dispersions were negligible relative to 
the galactic rotation velocity, then for disk lenses we would obtain a one-to-one relation 
v = 22Cte/(l — x) km/s , and thus the lens distance here would be D\ ens = 0.25L = 2.1 kpc. 
In reality, the random motions of both disk and bulge stars broaden this relationship 
somewhat, but we can still obtain a useful constraint. 

Given the observed v, we obtain a likelihood function 



L(x; v) oty/x{\- x) p L {x) v(l - xf / fsiys) - ^)( v © + v) + xw s ) dv s , (7) 



where pi is the density of lenses at distance x, and the integral is over combinations of 
source and lens velocity giving the observed v. and = (1 — x)(vq + v) +xv$ are the 
2-D source and lens distribution functions (normalized to unity). We adopt a disk velocity 
dispersion of 30 km/s in each direction, with a flat rotation curve of 220 km/s. We adopt a 
bulge velocity dispersion of 110 km/s, and no bulge rotation. For the density profiles, we 
use a standard double-exponential disk, and a barred bulge as in Han & Gould (1995a). 
The source is assumed to reside in the bulge, and the lens distribution functions are the 
sum of those for disk and bulge. 

The resulting likelihood function for -Di ens = X L is shown as the solid curve in Fig. 3; 
it is comparable to Figure 6 of Han & Gould (1995a), and is quite insensitive to specific 
parameter choices. This shows that the lens is much more likely to belong to the disk 
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than the bulge. The dashed curve is the likelihood function under the assumption of a 
main sequence lens. It is the direct result of multiplying the solid curve by the Gaussian 
probability that the lens brightness is consistent with the ± 4% upper limit, so its lower 
amplitude is an indication that much of the region allowed by the phase space constraint 
is not consistent with the main sequence lens brightness limit. 

It is interesting to note that inclusion of the constraint on the brightness of a main 
sequence lens changes the likelihood function significantly. The most probable distance 
from phase space and density considerations only is -Di en s = X L & 2 kpc, but this is excluded 
by the main sequence lens brightness constraint. There are two possible interpretations 
of this: either the lens is a main sequence star more than 2 kpc away, which is somewhat 
disfavored by the phase space constraint, or the lens is not a main sequence star. In the later 
case the lens could be an old white dwarf, or a neutron star, and it would be located at the 
distance indicated by the solid curve: D\ ens = 1.7^'^ kpc with a mass of M = 1.3+JjM 
(from eq. (6)). A significant population of stellar remnant lenses might also help to explain 
the number of long timescale lensing events that have been observed (Bennett et al. 1995, 
Han & Gould 1995b). If the lens is a main sequence star, we must use the dashed curve, 
which includes the lens brightness constraint and implies -Diens = %L = 2.8^'gkpc and 
M = O.49^Q23M . The relative normalization of the dashed curve in Fig. 3 indicates that 
a main sequence lens is somewhat improbable at the ~ 70% confidence level. 

An additional constraint we can place on the mass and location of the lens is due to 
the fact that, while the event has a high amplification of A max = 10.4 and the source has 
a relatively large radius of ~ 15 Rq, we cannot detect deviations from the point-source 
approximation. A light curve fit allowing for a finite source size indicates that the angular 
Einstein ring radius must be more than 10.2 times larger than the source. Using eq. (6), 
this implies that -Diens < 7 kpc, a result which is also implied by the likelihood function. 

If it were possible to make parallax measurements for a large fraction of the observed 
microlensing events, it would be possible to resolve most of the uncertainties in the in- 
terpretation of present microlensing results. Toward the LMC, for example, there appear 
to be more events than can be accounted for by known stellar populations (Alcock et al. 
1995a), but it is not known whether this is due to Machos in the halo, the disk, or in the 
LMC. Microlensing parallaxes would yield distance estimates, such as those presented here, 
which would be able to resolve this question for most events. Similarly, the microlensing 
optical depth toward the Galactic bulge seems to be higher than predicted (Udalski et 
al. 1994, Alcock et al. 1995b, Bennett et al. 1995). Suggested explanations of the high 
optical depth include a massive or "maximal" disk, a massive galactic bar pointed nearly 
along the line of sight (Paczyhski et al. 1994, Zhou, Rich & Spergel 1995), or a substantial 
fraction of the lensed sources residing behind the disk (Evans 1994). While the error bars 
on the position determinations from individual parallax events will not always allow us 
to distinguish between some of these different possibilities, parallax measurements for a 
substantial number of events will allow us to determine the location and mass function of 
the bulk of the lenses. 

Unfortunately, the prospects are poor that parallax measurements can be obtained for 
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many other microlensing events from the ground. The event described here is the longest 
of about 50 microlensing events that have been observed so far. In addition, its peak 
amplification occurs close to the middle of the bulge observing season. In the same data 
set (Bennett et al. 1995), we also have two events with t ~ 150 days but with the peak 
amplification occurring near the ends of our observing season, so any asymmetry is weakly 
constrained. We have also done parallax fits for these events, but obtain results with, at 
best, only marginal statistical significance even though the best fit projected transverse 
velocities are smaller than 50km/s. This situation could be improved somewhat in the 
future by obtaining more accurate photometry as events are discovered in progress (Alcock 
et al. 1994, 1995c; Udalski et al. 1995), but it is not possible to measure ground-based 
parallaxes for most of our events which have i values between 10 and 50 days? This means 
that ground-based parallax measurements will not allow us to resolve the unanswered 
questions about the properties of the bulk of the lensing objects. However, parallaxes for 
almost all events could be measured by a small satellite in solar orbit, as suggested by 
Gould (1994b, 1995). 

To summarize, we have presented the first detection of parallax in a gravitational 
microlensing event? We have used the fit parallax event parameters along with modest 
assumptions about the phase space distribution of lensing objects in a likelihood analysis 
to estimate the distance and mass of the lens. If the lens is near the distance indicated 
peak of the likelihood function, then it cannot be a main sequence star because we have 
not seen much light from the lens. We conclude that the lens is probably either a stellar 
remnant or located further away than about 2 kpc, which is consistent with, but slightly 
disfavored by, the likelihood analysis. Finally, we've considered the possibility of obtaining 
parallax events for the majority of events and concluded that this requires an additional 
telescope in space. 
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FIGURE CAPTIONS 



1. The observed two color light curve is shown as ±ler error bars, in linear units nor- 
malized to the fit unlensed brightness. The upper panel is the MACHO-R band data 
and the lower panel shows the MACHO-B band data. The dashed curve shows the 
constant velocity microlensing fit, and the solid curve shows the best-fit light curve 
allowing for the orbit of the Earth. Time is in days from JD 2449000. 

2. The lens mass is plotted as a function of the lens distance using the vt value deter- 
mined by the light curve fit shown in Fig. 1. 

3. Likelihood functions for the lens distance. The solid curve indicates the likelihood 
function incorporating the projected velocity and the Galactic model, while the 
dashed curve also includes the upper limit on the brightness of a main-sequence 
lens. 
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